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We study numerically the paramagnetic phase of the spin-1/2 random transverse-field Ising chain, 
using a mapping to non-interacting fermions. We extend our earlier work, Phys. Rev. 53, 8486 
(1996), to finite temperatures and to dynamical properties. Our results are consistent with the idea 
that there are "Griffiths-McCoy" singularities in the paramagnetic phase described by a continuously 
varying exponent z[S), where 5 measures the deviation from criticality. There are some discrepancies 
between the values of z(6) obtained from different quantities, but this may be due to corrections 
to scaling. The average on-site time dependent correlation function decays with a power law in the 
paramagnetic phase, namely r -1 / 2 '' 5 ', where r is imaginary time. However, the typical value decays 
with a stretched exponential behavior, exp( — ct 1 ^), where /i may be related to z(5). We also obtain 
results for the full probability distribution of time dependent correlation functions at different points 
in the paramagnetic phase. 



I. INTRODUCTION 

Most critical points occur at a finite temperature and 
one tunes through the transition by varying the tempera- 
ture itself. There are however, many critical points which 
occur at zero temperature and are traversed by varying 
some other parameter. These quantum critical points 
have recently been of great interest, particularly for sys- 
tems with disorder. One reason for this interest is that 
even the paramagnetic phase can have quite dramatic 
singularities. For example, in systems with a discrete 
symmetry, such as the Ising model discussed in this pa- 
per, there is a line in the phase diagram, on either side of 
the critical point, where correlations in time (but not in 
space) decay with a power lawETu. Power law behavior 
is characteristic of a critical point, but here this critical- 
like behavior occurs only in the time direction. We could 
say that there is a line of "semi-critical" points. Further- 
more, the exponent characterizing the power law. decay 
is expected to vary continuously along the lineOrti As a 
result static response functions may actv^.lb/j-diverge in a 
finite region away from the critical pointBotx 

These singularities arise from regions of the sample 
which have stronger interactions than average pand were 
first discussed, many years ago, by GriffithsO, in the 
context of classical models, where, however, they are 
rather weaku. At about the same time, McCoyd deter- 
mined exactly some properties of a two-dimensional clas- 
sical model (equivalent to the disordered one-dimensional 
quantum magnet that we study here), finding that the 
susceptibility diverges before the critical point is reached. 
It is now understood that this behavior is due the rare re- 
gions, more strongly coupled than average, discussed by 
Griffiths, but which give a larger effect in the quantum 
regime than near a classical transition. Hence we shall 
refer to Griffiths singularities in the quantum regime as 
Griffiths-McCoy singularities. 

We shall study here the disordered spin-1/2 Ising chain 



in one-dimension, for which many properties can be 
worked out in detail. The ground state of this model 
is closely related to the finite-temperature behavior of 
a two-dimensional classical Ising model with disorder 
perfectly correlated, along one direction, which was first 
studied by McCoyO and by McCoy and Wula. Subse- 
quently,|-the quantum model, was studied by Shankar and 
MurphyEI, and, in great detail, by D. S. FisherU From a 
real space renormalization group analysis, which becomes 
exact on large scales, Fisher obtained many new results 
and considerable physical insight. One of the key conclu- 
sions is that many quantities have very broad distribu- 
tions, with average and typical values being quite differ- 
ent. Consequently, a lot of information is lost by averag- 
ing. Confirmation of Fisher's surprising predictions for 
the T = 0, equal time behavior of the quantum problem, 
as well as some new results for distributions of various 
quantities, were obtained in earlier workEJ (henceforth 
denoted by YR) , which used a mapping of the spin prob- 
lem to non-interacting fermionsEjtJ to obtain accurate 
numerical results for quite large systems. 

In this paper we extend these techniques to finite tem- 
peratures and finite times. Our data suggest that the 
singularities are governed by a continuously varying ex- 
ponent, though there are some discrepancies which we 
discuss. We also obtain results for the distribution of 
time dependent correlation functions, showing, for ex- 
ample, that a typical correlation function decays with a 
stretched exponential dependence on time, as opposed to 
the average which decays with a power, as noted above. 
Recently, results for dynamical correlations at the criti- 
cal point have been obtained by Rieger and Igloitj using 
similar techniques to ours. Hence our results will be re- 
stricted to the paramagnetic phase. 



II. THE MODEL 

The model that we study is one-dimensional random 
transverse-field Ising chain with Hamiltonian 
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hiaf 



(1) 



Here the {erf} are Pauli spin matrices, and the interac- 
tions Ji and transverse fields hi are both independent 
random variables, with distributions 7r(J) and p(h) re- 
spectively. The lattice size is L, and, in this paper, we will 
take free, rather than the more usual periodic, boundary 
conditions. We will see later why this is necessary. 

Since the model is in one-dimension, we can perform 
a gauge transformation to make all the Ji and hi posi- 
tive. Unless otherwise stated, the numerical work used 
the following rectangular distribution: 



_. , j 1 for < J < 1 
^ ' ' otherwise 



p(h) 



h 1 for < h < h 
otherwise. 



(2) 



The model is therefore characterized by a single control 
parameter, ho- Defining 



A/j = [In/i] 
A. 7 = [lnJ] 



av 
av 



(3) 



where [. . .] av denote-an average over disorder, the critical 
point occurs wherota 



A h = A j 



(4) 



Clearly this is satisfied if the distributions of bonds and 
fields are equal, and the criticality of the model then 
follows from dualityu. A convenient measure of the devi- 
ation from criticality is given byu 



A, - A., 



var(ln h) + var(ln J) ' 



where var(- ■ •) denotes the variance. For the distribution 
in Eq. (0), we have 



1 



In ho 



(6) 



III. THE METHOD 

The numerical technique for static quantities, has been 
discussed in detail for periodic boundary conditions by 
YR. Here we consider free boundary conditions, which is 
simpler, so will only summarize— the main results. Fol- 
lowing Lieb, Schultz and Mattistil we start by using the 



the Jordan- Wigner transformation, which relates the spin 
operators to fermion creation and annihilation operators, 
cj and &: 



of =a\ + a% 

of = i(<4 ~ a i) 

of = 1 - 2a] ai = 1 - 2ck , 



(7) 



where 



4 = 4 ex P 



ai = exp 



i-i 



3=1 
i-l 



-VK c\c 



3 = 1 



(8) 



The Hamiltonian, Eq. (0), can then be writtenQ 

L L-l 



(9) 

Note that the fermion Hamiltonian, Eq (||), is bi- linear 
and so describes free fermions. 

We define operators ^\ for 1 < j < 2L by = cj 



and * 



t 

i+L 



where 1 < i < L. Similarly is the 



Hermitian conjugate of V E'^. Note that the ^>i and 
satisfy the usual fermion commutation relations for all i 
and j. 

The Hamiltonian, Eq (|9|) , can then be written in terms 
of a real-symmetric 2L x 2L matrix, H , as 



where H has the form 
H = 



A B 

-B -A 



(10) 



(11) 



(5) where A and B are L x L matrices with elements 



Ai.i 


= h L 


Ai.i+i 


= -Ji/2 


Ai+i.i 


= -Ji/2 


-Bi,i+1 


= Ji/2 


-Bi+l,i 


= -Ji/2 , 



(12) 

where, since we have free boundary conditions, elements 
with an index L + 1 are zero. Note that H is real sym- 
metric. 

Next we diagonalize H numerically, to find the single 
particle eigenstates with eigenvalues e M /2, /i = 1, 2, . . . 2L 
and eigenvectors ^ which are linear combinations of the 

with real coefficients0. It is easy to see that the eigen- 
states come in pairs, with eigenvectors that are Hermitian 



conjugates of each other and eigenvalues which are equal 
in magnitude and opposite in sign. We can therefore de- 
fine $t = -yt if Cfi > o and $L = 7^ if y! is the state 
with energy — e M . 

The Hamiltonian can then be written just in terms of 
L (rather than 2L) modes as 



L 



(13) 



where all the are now taken to be positive. The average 
energy per site is therefore given by 



1 L 



where n M is the Fermi function, 



1 



M exp(/3e M ) + 1 ' 
and the specific heat is consequently given by 

1 L 



(14) 



(15) 



(16) 



in units where fcs = 1, which we take from now on. 

We next consider equal time correlation functions, de- 
fined by 



(17) 



Since Sji = Sij, we can take j > i, without loss of gener- 
ality. Sij is given in terms of a determinant of size j — i 
byO 



Si 



where 



Gi,i+1 Gi,i+2 ' ' ' Gij 

Gj | ,.J • • • Gj-lJ 



G« = ((cJ-c i )(c] + c i )>. 



(18) 



(19) 



can be expressed in terms of the eigenvectors of 



the matrix H in Eq. (|ll[). Let us write 

L 
L 



where "0 and can be shown to be orthogonal matrices. 
Then 

Gij =((4-Ci)(ct+c,)) 



(21) 



(22) 



(23) 



At zero temperature, Eq. (|2l| ) goes over to Eq. (54) of 
YR. 

We now discuss how these results-apa generalized to 
time dependent correlation functions" 8 H 19 1. We are inter- 
ested in the a z -a z imaginary time correlation function, 
Sij(r) for < r < (3, where 



(lilt) = (7m7m> = 



(lllti) = 1 - (7 M 7u) = 



^•(r 1 -T 2 ) = ( < 7f(T 1 )a|(r 2 )), 



with 



Note that 



(24) 
(25) 
(26) 



which follows by cyclically permuting the trace, and so, 
without loss of generality, we just consider j > i. 
Substituting the transformation in Eq. (0) one has 



Sij(r) = ( exp 



m— 1 



x exp 



(ctW + Ci(r) 



i-i 



1=1 



This can be simplified since 

exp [-iirc^Cm] = A m B v 

where 



(27) 



(28) 



(29) 



and so 



i-l 



i-i 



(20) S y (r) = ( [[] A m (r)B m {T)]Mr)[l[ A^A, ) . (30) 



Hence Sij(r) involves the product of 2(i + j — 1) Fermi 
operators. This number is very large if one is interested 
in two sites in near the center of a large lattice, even if 
those sites are close together. 

The situation is much simpler for t — because then 
the product of all the operators to the left of site i is 
unityEJ, and so 

= /b<[ n ^iBi]aX (31) 

\ 2=i+l / 

where theJast line follows because A? = 1. As shown by 
Lieb et al.til, Wick's theorem, together with the observa- 
tion that 



(AiAj) = -(BiBj) =Si 



(32) 



enables one to write Sij(t — 0) as the Toeplitz determi- 
nant of order j — i in Eq. (O) . This is convenient because 
the determinant is small if i and j are close together, even 
for an sites far from the boundary of a large system. 

For r ^ in Eq. (^0|) one can still use Wick's the- 
orem but now there are many more pairs of operators 
to be included in the contractions. Wick's theorem for 
fermions requires the sum over all possible products of 
pair-averages, with a sign which is 1 or —1 depending on 
whether an even or odd permutation of the operators is 
necessary to get the operators in the product back to the 
original order. Thisjs called a Pfaffian, see e.g. the book 
by McCoy and Wv£j. If the number of operators, 2n, 
is large, (here n = i + j — 1), evaluation of the Pfaffian 
is intractable because the number of terms, (2n — 1)!!, 
is too large even for fast computers. However, the Pfaf- 
fian is also the square root of an antisymmetric matrix 
(of order 2n) formed from the pair averages. To be pre- 
cise, if A, B, ■ ■ ■ , Z are Fermi operators, and the average 
is over a free FermLHamiltonian (in the grand canonical 
ensemble), one hasE3 



(ABC ■■■Z) = 



(AB) (AC) 
-(AB) (BC) 

-(AC) -(BC) 

-(AZ) -(BZ) -(CZ) 



(AZ) 
(BZ) 
(CZ) 





1/2 



(33) 



As a simple example, consider i — 1, j = 2 in Eq. ( pQ ) 
for which Wick's theorem gives 

Si, 2 (r) = (A 1 {r)A 1 B 1 A 2 ) 

= (A 1 (t)A 1 )(B 1 A 2 ) - (A 1 (t)B 1 )(A 1 A 2 ) 

+ (A 1 (t)A 2 )(A 1 B 1 ). (34) 

This is easily shown to equal 



(Ax(r)Ax) (Ax(t)Bx) (A x (t)A 2 ) 1/2 

-(Ax(t)Ax) (AxBx) (AiA 2 ) 

-(Ai(t)Bi> -(AxBx) (B X A 2 ) 

-(A 1 (t)A 2 ) -(A,A 2 ) -(B,A 2 ) 

(35) 

which is just Eq. (^) for this case. Note that Eq. ( p2] ) 
gives (AiA 2 ) = and we also have A\ = 1. Hence, for 
t — 0, we get, from Eq. (j34|), the simpler result 



Sx,2(0) = (BxA 2 ), 



(36) 



which also follows immediately from Eq. JT§|). 

It is much more convenient to work evaluate the deter- 
minant and take the square root, than to directly evalu- 
ate the Pfaffian. In numerical work, the number of op- 
erations required to evaluate the determinant is of order 
(2n) , which, for n a few hundred, is feasible, whereas the 
(2n — 1)!! operations to evaluate the Pfaffian is definitely 
not. 

The pair averages needed for Wick's theorem can be 
evaluated in the same manner used to derive Eq. ( |2l| ) 
above for the equal time correlation functions. The result 
is, 



(Mtx)B {t 2 )) 
(B^Mti)) 
(Mtx)^^)) 



L 

ti=X 
L 



L 



/'./ ■ 



where 



V»(t) 



(1 



(37) 



(38) 



and t = 7~i — t 2 . From these pair averages the deter- 
minant in Eq. ( |33"| ) [with the operators in Eq. (|30|)] is 
evaluated numerically, and finally the square root taken. 
Since the imaginary time correlation function is real and 
positive, there is no ambiguity about the sign. 

We concentrate on the local correlation function, 
Su(t). This determines the local susceptibility from 



Xi 



Su(t) dr. 



(39) 



We determine Su(t) for different values of t, increasing 
in a roughly logarithmic manner, and approximate the 
integral by the trapezium rule. We compute the average 



local correlation function and average local susceptibility, 
defined by 



Sav(r) = 7$>«( T ) 



= 7 (40) 



From now on, for compactness of notation, we will de- 
note Su(t) by S(t). In addition, because there are large 
fluctuations in the values of S(t) from site to site, we also 
look at the distributions of this quantity for different r. 



IV. PHENOMENOLOGICAL DESCRIPTION 

From recent work,0~B0El a phenomenological descrip- 
tion of the Griffiths-McCoy region of the paramagnetic 
phase has emerged. Singularities arise from small regions 
which are "locally in the ferromagnetic phase" and have 
a very small energy gap. As a result, there are low en- 
ergy "cluster" excitations, which have a power law distri- 
bution of energies, AE. The probability of having a low 
energy excitation is proportional to the size, L, and so we 
can write the dimensionless probability, AE P(AE),a,s 



AEP(AE) ~ L(AE) 



l/z(S) 



(41) 



where we write the power in terms of a dynamical ex- 
ponent z(8), since Eq. ( pT| ) corresponds to the standard 
relation between a length scale (L here) and an energy 
scale (AE) here. The exponent is expected to vary con- 
tinuously in the Griffiths-McCoy phase, and we indicate 
this by the notation, z(8). 

For the distribution used here, Eq. (||), in the limit 
8 — ► oo we can neglect the interactions and the excitation 
energy of a single spin is just 2hi. Since the hi have a 
uniform distribution for h — ► 0, it follows that z = 1 in 
this limit, i.e. 



lim z(8) = 1. 

S — >oo 



(42) 



Furthermore, it has been establishedi'i that z(8) oo 
at the critical point, i.e. 



lim z{8) = oo. 

5-.0 



(43) 



The energy per site at low temperature can be esti- 
mated from the energy of the excited clusters, i.e. 

EAT) - Kv(0) « \ [ AEP(AE) -^tf^L 
L J 1 + exp (— pAE) 



^ rpl+l/z(S) 

so the specific heat varies as 
for T -» 0. 



(44) 
(45) 



For a site in a cluster with a low energy excitation AE, 
the long time behavior of S(t) is exp(— AEt). Averaging 
over the distribution gives 



Sb, v (t) 



1 



T l/z{5) ' 

and integrating r up to (3 gives 

j>-l+l/z(r5) 



loc 

Aav 



(46) 



(47) 



for T -> 0. 

To our knowledge, not much is known about the dis- 
tribution of S (t) and one of the goals of the present work 
is to deduce its form. Note that because we expect the 
distribution to be very broad, with average and typical 
values quite different, we cannot predict its form simply 
from knowing the average. 



V. RESULTS 
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FIG. 1. The specific heat as a function of T for 
8 = 0.549 (fto = 3) and sizes L = 64 and 32. According 
to Eq. @, the slope is l/z(5). 

We start with our results for the specific heat. The 
data for ho = 3 (8 = 0.549), shown in Fig. I], has good 
straight line behavior for more than five decades, the 
slope, equal to l/z(8), is 0.74. 

According to Eq. p2] ) z(8) should tend to unity as 
8 — > oo. This is confirmed by the data for the specific 
heat for h = 50 (8 = 1.956), shown in Fig. |2[ which has 
a slope of 0.99, very close to the expected value. 

The temperature dependence of the susceptibility is 
shown in Fig. || for ho = 2 (8 = 0.346). Again a power 
law behavior is obtained, as expected from Eq. (l47|). 
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FIG. 2. The specific heat as a function of T for 
8 — 1.956 (ho = 50) and sizes L = 64 and 32. According 
to Eqs. ( E^ ) and (j42J) , the slope should be close to unity, as 
indeed it is. 



We now proceed to our results for the time dependent 
correlation functions. Data for 5 av (r) for h = 2 and 
6 (5 = 0.346 and 0.895) at T = are shown in Figs. | 
and H on a double logarithmic scale. The results clearly 
indicate a power law behavior, as expected from Eq. (|4(f). 
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FIG. 4. The average on-site (imaginary) time dependent 
correlation function at T=0 for 8 = 0.346. The lattice size is 
L = 64. 
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FIG. 3. The susceptibility heat as a function of T for 
8 = 0.346 and size L = 64. 



FIG. 5. As for Fig. |j, but with 8 = 0.895. 

From our results for C av and S av (r), jsse obtain z(S) 
from Eqs. (45) and (^6|). We summarizes these results 
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FIG. 6. The dynamic exponent z(S) obtained both from 
the specific heat and the time dependent correlation func- 
tion. The short dashed line indicates the asymptotic value 
z{5 — >rpo) = 1, and the long dashed line is the prediction of 
Fishera, z -1 = 28 for 5^0. At the critical point, 5 = 0, z{5) 
is predicted to be infinite. 
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FIG. 7. The distribution of - In S(t) for 5 = 0.895. The 
lattice size is L = 64 and 50, 000 samples were averaged over. 



in Fig. ||. While the trend in the two sets of data is the 
same, there are some differences, which we do not un- 
derstand very well. The lattice sizes used are quite large 
and the data used to generate the estimates for z(6) fit 
a straight line over a fairly large range, especially for the 
specific heat. This suggests that the estimates for z{§) 
should agree well. However, corrections to scaling may 
be large because free boundary conditions are used here, 
as opposed to the more conventional periodic boundary 
conditions. 

So far we have just considered the average value of 
various quantities. However, one of the most surprising 
features of this model is that distributions are so broad 
that average and typical values can be quite different. We 
have therefore also studied the distribution of the on-site 
time dependent correlation functions for different values 
of r. 

Fig. ^ shows the distribution of lnS(r) for 5 = 0.895. 
The distribution is broad and becomes broader with in- 
creasing r. We therefore attempt to collapse the data by 
plotting it versus the scaling variable 



In S(t) 
T x/n 



(48) 



with an appropriate choice for \i. We shall see that, at 
a given value of 5, the distributions of S(t) for different 
times are all given by a single scaling function P(x). 

Fig. H shows that the data plotted in this way collapses 
quite well with the choice /i = 2.2. For small x the behav- 
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FIG. 8. A scaling plot of the data in Fig. [?]. The scaling 
variable is given in Eq. ( |48[ ) and we have taken ^ — 2.2 to get 
the best fit. The solid line is Eq. ( |5C| ) with z = 1.2, consistent 
with Eq. @, and c = 0.2. The fit to the data is good. 



ior seems to be close to linear, though the data for small 
times has a finite intercept, which perhaps vanishes for 
t — > oo. We shall discuss this again later. Scaling plots 
for S = 0.549 and 0.346 are shown in Figs. || and The 
values of fj, used in these fits are 2.2 and 2.5. 

From a variable range hopping analysis, 
M. P. A. FisherBfl has argued that \x should be related 
to z(5) by 



l + z(S). 



(49) 



In Appendix A we show that, within the same set of 
assumptions, the scaling function is given by 



P{x) = c(cx) x l z exp 



[ex) 



(z+l)/z 



where 



- In S(t) 

-l/(z+l) ' 



(50) 



(51) 



as also follows from Eqs. (|4^) and (|4^), and c is a scale 
factor. The data for S = 0.895 fits the scaling function 
in Eq. ( |50| ) quite well, with an appropriate choice of the 
scale factor c and a value for z[5), and hence /i, consis- 
tent with that from other data shown in Fig. |6[ How- 
ever, the data closer to the critical point in Figs. || and 
|lO| only fit the scaling form if z[5) in Eq. ( p30| ) and \x in 
Eq. ( [48|) are adjusted independently. The values for z[5) 
are found to be in reasonable agreement with those in 
Fig. |6| but the values for /i are then somewhat incon- 
sistent with Eq. (|4S|). NoteJaowever, that the stretched 
exponential is only expectecO to be valid for times larger 
than a characteristic time which diverges at the transi- 
tion. Hence it is possible that our data is not at suffi- 
ciently long times for smaller S to get a good estimate for 
M- 

The fact that the scaling variable is given by Eq. (|4g| ) 
shows that the typical correlation function, which we de- 
fine to be the exponential of the average of the log, is 
given by 



S typ [r) = exp([lnS(r)] av ) ~ exp(-cr 1 ^) 



(52) 



where c is a constant, i.e. a stretched exponential. This 
contrasts with the average, which varies as a power of r, 
as shown in Eq. (f46|). 

The average value is dominated by rare regions which 
have a much larger correlation function at long times 
than in a typical region. It is interesting to ask whether 
the average value is contained within the scaling function, 
or whether it comes from contributions which are actually 
corrections to scaling. We shall argue that there are both 
scaling and and non-scaling contributions to the r 1 / 2 ^ 5 ) 
behavior in Eq. (f46|). Similar behavior has been found 
recently for the distribution of the equaktime end-to-end 
correlation function at the critical pointed. 

To estimate the contribution from the scaling function, 
we assume that the distribution, P[x), of the scaling vari- 
able x, in Eq. (^8|), has the form x a in the limit x — » 0, the 
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FIG. 9. As for Fig. | but with 5 = 0.549. The best fit 
is for fi — 2.2. The solid line corresponds to the scaling dis- 
tribution in Eq. ( [so]) with z(8) — 1.36 and the scale factor 
c = 0.36. Although the fit is not too bad, it is actually some- 
what inconsistent because z{5) and /i should be related by 
Eq. (^i|), which would require z — 1.2. However, this choice 
of z(S) works less well. 
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FIG. 10. As for Fig. g but with 8 = 0.346. The best 
fit is for /i = 2.5. The solid line corresponds to the scaling 
distribution in Eq. ( jHoj ) with z{8) — 2.0 and the scale factor 
c = 0.46. However, this is inconsistent since, according to 
the theory in the appendix, n and z(S) should be related by 
Eq. (^9|), which would require z = 1.5. However, this choice 
of z(S) works less well. 



only region which contributes to the average. This is in- 
deed the case for the phenomenological theory discussed 
in the appendix, for which P(x) is given in Eq. (50). 
Since S(t) — exp(— it 1 /''), it follows that the scaling 
contribution to the average is 



/>OC 

Jo 



1 



r(l+0)//i' 



(53) 



which is of the form in Eq. (Eq) with the identification 



z(S) 



/' 



1 + a 



(54) 



Note that the phenomenological theory, Eqs. (^)-(^Tj), 
satisfies this condition with 



1 



z(S) 



(55) 



Although the numerical data only fits the phenomeno- 
logical theory close to the critical point with values of /j, 
and z(S) which are inconsistent with Eq. ((49|), this may 
be due to the times studied not being long enough and 
other corrections to scaling. It seems very likely that 
there is a contribution to the average time dependent 
correlation function from the scaling function. 

However, we shall now argue that there is an addi- 
tional contribution to the average correlation function, 
which comes from corrections to scaling. The exponent 
a is positive, so P(0) = in the scaling limit. How- 
ever, we see from the data in Fig. ^ that the value of 
Prob (— ln5(r) = 0) is finite for finite r though it does 
decrease for large r. This finite intercept is a correction 
to scaling. It is easy to see that if Prob (— In S(t) = 
0) ~ 7- -1 /*(*), then this finite intercept will give a power 
law contribution to the average correlation function of the 
form in Eq. @. Results for Prob (- In S(t) = 0) for dif- 
ferent values of r at 5 = 0.549 are shown in Fig. O. The 
slope is —0.55, which is a little less than the estimate for 
l/z(S) obtained by other methods, see Fig. 0, but since 
there are errors in extrapolating the data to S(t) = 0, 
this is perhaps not significant. Hence our results suggest 
that Prob (— ln5(r) = 0) vanishes as t~ x / z ( s \ which 
gives an additional contribution to the average correla- 
tion function which is not contained within the scaling 
function. 



VI. CONCLUSIONS 

We have studied the paramagnetic phase of the disor- 
dered one-dimensional Ising chain in a transverse field. 
Average properties (equal time, time dependent, and 
temperature dependent) can be characterized by a con- 
tinuously varying exponent z(S). At criticality, z(0) = oo 
and, for the model studied, lim^oo z{5) = 1. As an ex- 
ample, the average, on-site, time-dependent correlation 
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FIG. 11. The probability of y = -ln5(r) at y = for 
different values of r for 8 — 0.549 and L = 64. The slope 
is —0.55 which is somewhat less than, but not too far from 
the value of l/z(6) in Fig. |(| Hence this finite value of the 
probability at lnS(r) = (which is a correction to the scal- 
ing form) gives a contribution to the power law decay of the 
average correlation function shown in Eq. 



function decays with a continuously varying exponent, 
see Eq. ([46|). There are some discrepancies between the 
values of z(S) obtained in different ways. We suspect that 
they are due to corrections to scaling, but one should per- 
haps also worry that the simple phenomenological picture 
in Sec. IV may be inadequate. 

By contrast, the typical time dependent correlation 
function decays with the stretched exponential form in 
Eq. (p3), where the exponent fi is probably given by 
Eq. (|49|), though our data close to the critical point shows 
some discrepancy with this. At long times, the distribu- 
tion of S(t) has a scaling form, being a function, P(x), 
of the single variable x in Eq. pq). This may be given 
by Eq. (p0|), though corrections to scaling in the small 
x region prevent us from checking this thoroughly. For 
example, the data in the scaling plots in Figs. p[]lO| are 
expected to go through the origin, but they actually do 
not, the discrepancy being larger for smaller r. 

The average of S(t) arises both from the small x region 
of the scaling function, and from corrections to scaling. 

It is interesting to ask to what extent the behavior in 
the paramegnetic phase is universal. Presumably the de- 
pendence of p. on 5 is non-universal (except close to the 
critical points) but, for a given z, is the scaling function 
for the ditribution of local correlation functions univer- 
sal? Since the scaling function involves the limit of long 
times, it is possible that the microscopic details do not 
matter, only the form of the low energy density of states, 



in which case the distribution would be universal. 

It would be interesting to see to what extent the results 
found here in one dimension go over to higher dimensions. 
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APPENDIX A: 

Following the suggestion of M. P. A. Fishcr,0 we as- 
sume a simple phenomenological picture in which there 
are cluster excitations localized about different sites . It 
is plausible that the wavefuntion of an excitation will de- 
cay exponentially with a distance of order the typical cor- 
relation function, £typ- The contribution to the dynamics 
of site at the origin from an excitation centered n lattice 
spacings away will therefore be exp (— r/£ t y P — £«■?")■ The 
local correlation function is therefore given by 



5(r)=A^exp 



where 



n=0 



A = 1 - e 1 ^ 



r 
£typ 



(~ r 1 for e » i) 



(Al) 



(A2) 



ensures that S(0) = 1. 

As a simple model we will assume that the e„ are un- 
corrected. Writing e n = eoe„, where eo is a characteristic 
energy scale, we take the distribution of the e n to be 



p(e) = (A + l)e\ (0<e<l), 
where we see from Eq, ( ffl| ) that 

A = 1 ' 



(A3) 



(A4) 



Note that at large r there i s a c ompetition between the 
exp(— r/^typ) factor in Eq. (Al), which decreases with 
increasing r and so would prefer to have r small, and the 
exp(— e„r) factor which wants to have the smallest e n , 
which may be on a site far away. 

The aver age value is easy to work out, since each term 
in Eq. (Al) can be evaluated separately, with the result 



[S(r)] av = 



r(i + i/z) 



(eor)V 

for t» 1, where we used Eq. (|A 



(A5) 



However, the average is very different from the typical 
behavior. To see this, we will determine the full distribu- 
tion of S(t). The major simplification in the calc ulation 
is that, at large times, the exponent in Eq. ( |Al[ ) varies 
over a large range, so the sum wi ll b e dominated by the 
largest single term. Writing Eq. ( |Al[ ) as 



S(t) 



A^exp 

n=0 



Xn 
styp 



where 

x n = n + $typ£oTe„, 
then S(t) is given for large r by 

S(t) « A exp 



£ty P 



(A6) 



(A7) 



(A8) 



where x m i n is the smallest of the x n . 

For the time being we will work in units of time where 
£typ£o = 1. At the end, we will put back this factor by 
replacing r by CtypfoT- 

The probability of x n , is given by 



*Kn{Xn) 



A + l 



(A9) 



for n < x n < n + r, and zero otherwise. Notice that the 
different x n have different distributions. 

We now wish to determine the probability that mini- 
mum value of the x n is x m i n . Let us assume first that it 
is x m which has the smallest value. The probability that 
(i) x m is the smallest, and (ii) its value is x m i n , is given 

by 



where 



J^n (*^min) 
Qn (^min) 



i (^min) 



n+r 



7T n (x) dx. 



(A10) 



(All) 



Since ir n (x) vanishes for x < n, it follows that for 
n > x min , the full range over which n n (x) is non-zero 
is included in the last integral, and so 



t(-^min) — 1; (^ ^* *^min)- 



(A12) 



As a result, the product in Eq. ( A10 ) need only be taken 
up to the last integer below x m i n . For n < x m i n , one has 



i (Xn 



A+l 



(A13) 



We will see that x m i n <C r and the important values 
of to are those where 1 < m C r, and so we can rewrite 



Eq. (A10) as 




x exp 



(A14) 



Expanding the log using a; m i n — x -C t and performing 
the integral gives 



A + 1 fx 



exp 



P A+2 



(A + 2)r 



A+l 



(A15) 



Remember that this is the probability that x m is the 
smallest of the x n and its value if x min . We therefore next 
sum over all m less than x min to get the total probability 
that the minimum is x m i n , i.e. 



i(^min) 



A+l 



exp 



„A+2 



(A16) 



(A + 2)r A +! 

Replacing r by eo£typT, and chan ging variables from 



to y : 



ln(5(r)/A) using Eq. (AE), gives 



p y(y) = ftyp 



A+i 



exp 



A+2 



A + l(e r)^+i 



(A17) 



It is easy to check that this yields the average value in 
Eq. (|aH). 

Eq. (A17) can be cast in a scaling form if we define 



,1/A 



ln(S(r)M) 

r l/A 



(A18) 



Note, however, that the scaling limit involves taking 
— In5(r) — > oo so the lnA term in Eq. (A18) repre- 
sents an additive correction to scaling and can be omit- 
ted. This then leads to Eq. (pl|), in which A has been 
replaced by z using Eq, (A4). The probability of x is 
then a function just of x (apart from a scale factor) i.e. 



P{x) = c(cx) 1/z exp 



z + 1 



(car) 



(z + l)/z 



(A19) 



where we have again ex pres sed the result in terms of z 
rather than A using Eq. (A4), and 



(£typ e O 



Eq. (|A2q ) is precisely Eq. (|5C|) of the text 
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